Optical mesh lattices with PT- symmetry 



t-H 

o 

< 
00 



Mohammad- Ali Miri 1 '*, Alois Regensburger 2 ' 3 , Ulf Peschel 2 , and Demetrios N. Christodoulides 1 
1 CRE0L/ College of Optics, University of Central Florida, Orlando, Florida 32816, USA 
2 Institute of Optics, Information and Photonics, 
University of Erlang en- Nuernberg, 91058 Erlangen, Germany 
3 Max Planck Institute for the Science of Light, 91058 Erlangen, Germany 
(Dated: August 9, 2012) 

We investigate a new class of optical mesh periodic structures that are discretized in both the 
transverse and longitudinal directions. These networks are composed of waveguide arrays that 
are discretely coupled while phase elements are also inserted to discretely control their effective 
potentials and can be realized both in the temporal and the spatial domain. Their band structure 
and impulse response is studied in both the passive and parity-time (VT) symmetric regime. The 
possibility of band merging and the emergence of exceptional points along with the associated 
optical dynamics are considered in detail both above and below the PT-symmetry breaking point. 
Finally unidirectional invisibility in PT- synthetic mesh lattices is also examined along with possible 
super luminal light transport dynamics. 

PACS numbers: 42.25.Bs, 11.30.Er, 42.82.Et 
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I. INTRODUCTION 

Optical wave propagation in periodic structures has 
been a theme of considerable attention in the last twenty 
years or so [1-5]. In general such arrangements can 
exhibit intriguing and potentially useful light dynamics 
that are otherwise impossible in the bulk. Periodic 
Bragg gratings and two and three dimensional photonic 
crystals with a complete band gap are examples of such 
configurations. Optical waveguide arrays represent yet 
another class of periodic structures and have been the 
focus of intense research in the last decade [(Sl-Hi}. As 
indicated in several studies, optical arrays can be used 
as versatile platforms to observe a number of processes 
ranging from Bloch oscillations to Landau-Zener 

tunneling and from Anderson localization [|2| to 

discrete solitons 0, B a nd Rabi oscillations |l3| and 
dynamic localization [14[, just to mention a few. Along 
similar lines, longitudinally periodically modulated 
optical waveguide arrays have also been studied in 
several works [HI, Qij . 

Quite recently, the concept of parity-time (VT) 
symmetry has been introduced in the field of optics 
jl7H19|. Interestingly, the very idea of PT-symmetry 
originated within the framework of quantum mechanics 
through which a large class of complex non-Hermitian 
Hamiltonians was identified that could in principle 
exhibit entirely real spectra (2(| This can happen 
as long as the associated Hamiltonian and the combined 
VT operator share the same set of eigenfunctions. In 
this case, this is possible provided that the corresponding 
complex potential satisfies the condition V*(x) = V(—x) 
[20l | . This directly implies that the real and imaginary 
parts of the potential must be even and odd functions of 
position respectively. 

Lately, optical PT-symmetry has been experimentally 



observed in two-element coupled systems where non- 
reciprocal dynamics (disrupting left-right symmetry) 
and spontaneous PT-symmetry breaking has been 
demonstrated for the first time [H, [23|. What made 
this transition to photonics possible is the isomorphism 
between the evolution equations in quantum mechanics 
and optics. Indeed one can show that the complex 
refractive index, n(x) = tir(x) + m/(x), plays in this 
case the role of an optical potential. In this repre- 
sentation the real part ur(x) stands for the refractive 
index profile while the imaginary part ni(x) represents 
the gain or loss in the system (depending on its sign). 
Clearly, under PT-symmetric conditions one expects 
that tir(x) = ur(—x) and ni(x) = —nj(—x). In other 
words the index distribution must be an even function of 
position whereas the gain/loss must be anti-symmetric. 
Thus far, several works have pointed out that VT- 
symmetry can lead to altogether new optical dynamics 
which are otherwise impossible in standard passive 
optical arrangements [2J-[4l|. These may include for 
example the occurrence of abrupt phase transitions along 
with the appearance of the so-called exceptional points 
[24]-[26j, power oscillations [18j |. breaking left-right sym- 
metry and the occurrence of secondary emissions jl8| . 
In addition new classes of optical solitons pjl [38j and 
nonlinear VT optical isolators [29| , have been suggested 
along with unidirectional invisibility [35|, [36|, broad 
area VT single mode lasers [30], and coherent perfect 
absorbers [32l-l34j . Note that phase transitions similar 
to those occurring due to T^T-symmetry breaking have 
also been reported in other systems involving gain/loss 
modulation 14211 . 



So far however, experimental observations of VT- 
symmetry in optics have been carried out in basic 
coupled systems involving only two elements [H, |23| . 
What has hindered progress along these lines is not 
only the delicate balance needed between gain and 
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loss but also the requirement that the real part of the 
potential should remain symmetric, even in the presence 
of gain and loss. Fulfilling all these conditions at the 
same time is by itself a challenging task because of 
the underlying Kramers-Kronig relations. Therefore 
of interest would be to develop a new class of optical 
platforms where refractive effects and gain or loss can 
be treated separately and thus facilitate the realization 
of T^T-symmetric optics on a large scale. This goal is 
reached by discretizing not only the transverse, but also 
the longitudinal or propagation direction. 

In this work we introduce a new class of PT-symmetric 
optical lattices. These mesh arrangements are composed 
of an array of waveguides with each one of them being 
discretely and periodically coupled to its adjacent 
neighbors (Fig. Ufa)). Unlike ordinary waveguide arrays, 
light propagation in such mesh systems is discretized 
in two dimensions (transverse and longitudinal). The 
band structure of this family of mesh lattices is derived 
analytically and its effects on light dynamics are investi- 
gated. Because of the aforementioned 2D-discretization, 
the resulting band structure is characterized by both a 
transverse and a longitudinal Bloch momentum. 

As we will see, this type of lattice can provide a 
versatile platform for observing a host of T^T-symmetric 
phenomena and processes. Along these lines, phase 
elements can be readily inserted in the mesh lattice so 
as to control the real part of the array potential while 
amplifiers (that are turned on or off) can be included to 
provide the needed anti-symmetric gain/loss profile. The 
fundamental building block of such a mesh structure 
happens to be a basic PT-symmetric coupler arrange- 
ment. What makes this structure practically appealing 
is the physical separation between the coupling and 
amplification/attenuation stages within this building 
block. Band merging effects as well as the emergence of 
exceptional points are investigated in this family of VT 
lattices along with superluminal light transport. Finally, 
unidirectional invisibility in T^T- synthetic mesh lattices 
is also examined and pertinent examples are provided. 



II. OPTICAL MESH LATTICES IN THE TIME 
DOMAIN 

Lately, the temporal equivalent of an optical mesh 
lattice has been experimentally realized using time- 
multiplexed loop arrangements [43j . Such configurations 
have been systematically employed to investigate a 
number of issues ranging from discrete quantum walks 
[4i-[47j to Bloch oscillations and fractal patterns (HE^. 
While spatial realizations of such mesh lattices have 
also been reported [13, |49| , time- multiplexed fiber loop 
schemes have so far demonstrated a high degree of flex- 
ibility [43|, |45|. In time- multiplexed schemes, a discrete 



time axis n corresponds to the transverse discrete axis 
of a corresponding spatial optical mesh lattice as shown 
in Figs. DJc,d). These time- multiplexed configurations 
involve two coupled fiber loops-coupled via a central 
50/50 directional coupler. These two fiber loops differ 
in length by AL. Here, the equivalent transverse 
coupling to the left and right sites is enabled by this 
length difference between two loops. An independent 
discretization in time is then obtained by monitoring 
the round-trip number m in these loops. Hence, the 
system is essentially discretized in two-dimensions. As 
we will see, the propagation dynamics of light pulses in 
such discrete temporal lattices, are exactly identical to 
those expected in the spatial configurations discussed in 
section III of this paper. 




FIG. 1. (Color online) (a) Discrete mesh lattice in the spatial 
domain; (b) pulses in the spatial lattice propagate to the left 
and right on a discrete ID position grid; (c) equivalent time- 
multiplexed scheme involving two coupled fiber loops with a 
length difference AL; (d) pulses in the short (long) loop are 
advanced (delayed) in time, thus making discrete steps on a 
ID temporal grid. The propagation dynamics are identical to 
those expected in (b). 



In an experimental set-up, gain and loss can be 
readily integrated into the fiber loops, e.g. by standard 
semiconductor or fiber optical amplifiers and amplitude 
modulators. This in turn may enable the demonstra- 
tion of large-scale PT-synthetic optical lattices in the 
temporal domain [Hoj| . Given that the "topological" 
arrangements of a temporal and a spatial mesh lattice 
are totally equivalent, here, without any loss of general- 
ity we will consider for simplicity their spatial realization. 
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III. OPTICAL MESH LATTICE AND ITS BAND 
STRUCTURE 



Figure [2] illustrates the spatial realization of such 
a mesh lattice when only passive phase elements are 
involved. As previously indicated, this configuration 
can be synthesized using an array of waveguides that 
are periodically and discretely coupled to their next 
neighbors (at the rectangular regions of Fig. [2j). In ad- 
dition phase elements can also be inserted. Each phase 
element introduces at every array site n a phase (j) n that 
happens to be independent of the discrete propagation 
step m. The location of each phase modulator in the 
lattice is denoted in the figure by a circle. As we will 
later demonstrate, these phase modulators effectively 
play the role of a refractive index profile in continuous 
arrangements. Figure QJa) schematically shows how light 
flows in such a system when only one of the waveguides 
is initially excited. After traveling a certain distance 
in each waveguide, light couples to the adjacent left 
(right) channel through a coupler, and after propagating 
this same distance it then couples to the adjacent 
waveguide to its right (left). Indeed light propaga- 
tion in this system leads to an interference process that 
is equivalent to a discrete time quantum walk [3, HH, [52| . 




FIG. 2. (Color online) An optical mesh lattice; The lattice 
is composed of an array of waveguides which are periodically 
coupled together in discrete intervals. Circles indicate the po- 
sition of phase elements and rectangles the coupling regions. 
The dashed lines show the discrete points where the field in- 
tensity is evaluated before coupling occurs. 



As Fig. [2] clearly indicates, this mesh lattice is di- 
atomic in nature. Using the simple input /output rela- 
tion of a 50:50 coupler [53j and by considering the effect 
of the phase elements, it is straightforward to show that 
the light evolution equation in this system takes the form: 



p ic f>n 

m+1 c r / tyx i ■rm\ 

a n = — [K + «»n ) 



(la) 



■ —%0n ( rn i -7 m \i 

+ e (-a n -i +ib n _ 1 )\, 



e ict>n 

7 m+1 ° r / 7 ttl i -?i 



(lb) 



+ C*»+l(-C + l+KC +1 ) 



In Eqs. (pQ), a™ and b™ represent the field amplitudes 
at adjacent waveguide sites n (in the n'th column) at 
a discrete propagation step or distance m ( m'th row). 
It should be noted that in deriving these equations the 
phase accumulated due to propagation in any waveguide 
section is ignored. Indeed a waveguide section of length 
/ between two subsequent couplers leads to a phase 
accumulation of /3Z, where f3 is the propagation constant 
of the guide. Yet, one can readily show that even in the 
presence of these additional phase terms Eqs. (pp) remain 
the same once a simple gauge transformation is used; 
(a™ 6-)^(a™ ^)e j2m/3i . 

To establish the necessary periodicity, we assume that 
the phase elements provide a phase potential that alter- 
nates between two different values in n: 



+0o, n : even 
—0o, n : odd 



(2) 



This kind of phase potential has a translational symme- 
try n +2 = 0n which leads to a transverse periodicity 
in this "four-atom" lattice with a fundamental period 
of TV = 2 where each cell is diatomic. In addition the 
lattice is now periodic in both n and m. 

First we study the band structure of this mesh system. 
Once the band characteristics and corresponding Bloch 
modes are known, the dynamic properties of the system 
can then be extrapolated. To find the dispersion relation 
of this lattice we consider discrete " plane wave solutions" 
of the form e l Q n e ldm where Q represents a Bloch momen- 
tum in the transverse direction and plays the role of a 
propagation constant. To obtain the corresponding band 
structure we assume solutions of the form: 



B r 



(3) 



where A n and B n are periodic Bloch functions with the 
period of N = 2, i.e. A n+ 2 = A n and B n +2 — B n . In 
general, for n = 2j, we use A ni B n = Aq,Bq while for 
n = 2 j + 1 we employ A ni B n = Ai,B±. This comes 
from the fact that a unit cell of this periodic structure 
includes two discrete positions n. 



By inserting Eqs. (j3j) in (pQ), and by adopting the phase 
potential of Eq. (|2]), we obtain the following dispersion 
relation after expanding the corresponding 4x4 deter- 
minant of a unit cell: 

cos(2Q) = 8cos 2 (6)-8cos((j)o)cos(6)+4cos 2 ((i)o)-3 (4) 

As expected from the double periodicity of this system 
in both n and m the band structure is also periodic 
in both Q and 6 having fundamental periods of tt and 
2tt respectively. This represents a major departure from 
optical waveguide arrays where the propagation dimen- 
sion is a continuous variable. Under the assumption of 
Eq. (j2]), this mesh arrangement exhibits four primary 
bands which are periodic with respect to the two Bloch 
momenta. Figure [3] depicts the band structure of this 
mesh lattice when 6r, = 0.27T. 




FIG. 3. Band structure of the optical mesh lattice in the pres- 
ence of periodic step-like potential created from phases, alter- 
nating between — 0o and +0o where 0o = 0.2tv. The shaded 
area shows the band gap regions and the dotted boundary 
depicts the primary Brillouin zone of this lattice. 

Equation |4] is valid in general for any arbitrary choice 
of cV However it should be noticed that in the special 
case where 0o = (empty lattice) this relation becomes 
degenerate. Indeed for the empty lattice the periodic- 
ity of this diatomic lattice is TV = 1 and hence its Bril- 
louin zone involves two bands and lies in the domain of 
— 7r < Q < 7r and — tt < < tt. The folded version of 
this Brillouin zone (corresponding to the empty lattice) 
is shown in Fig. HJa) where the two bands are degener- 
ately folded into four. Figures |4jb,c,d) depict the band 
structure of this mesh lattice for three nonzero values of 
0o within the Brillouin zone as a function of the Bloch 
momenta, i.e., — tt/2 < Q < tt/2 and — tt < < tt. Again 
the shaded areas show the associated band gaps. Accord- 
ing to Fig. 3J a nonzero 0o lifts the degeneracy and leads 
indeed to four bands. 

According to Eq. ((4]) and as one can see from the figures 
the band structure has a reflection symmetry around Q = 
and = 0. For any finite 0o there are four bands in 
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FIG. 4. (Color online) Band structure of an optical mesh 
lattice for several cases; (a) Lattice without any phase po- 
tential 0o = (empty lattice), (b) Lattice with a symmet- 
ric phase step- like potential varying between — 0o and +0o 
when 0o = 0.27T, (c) same as in (b) but with 0o = 0.57T, (d) 
0o = 0.77T. For case (a) the reduced Brillouin zone is de- 
picted while for the rest the first Brillouin zone is shown in 
its entirety. 



the Brillouin zone, all having a zero slope at the center 
(Q = 0) and at the edges (Q = tt/2). For the empty 
lattice on the other hand, in reality there are two bands 
and the slope is zero at the center (Q = 0) of the top 
band while it is non-zero at the two edges (Q = tt/2) 
and at Q = 6 = where the bands collide and there is 
no band gap between them. The addition of the phase 
potential 0o to the empty lattice breaks this degeneracy 
and creates band gaps at these points. This breaking 
of the degeneracy becomes clear by comparing Figs. Ufa) 
and (b). Equation [4] can also be written in a more explicit 
form as a function of Q: 



±cos 



COSyCpO 



± 



\Jcos 2 (Q) + sin 2 



(5) 



where in this relation any combination of the two 
plus/minus signs corresponds to each of the four bands. 

Before ending this discussion, it is worth noting that 
this phase potential does not need to be symmetrized in 
a 0o fashion as done before in this section. In fact any 
periodic potential that is alternating in n between two 
different phase values will break the degeneracy of an 
empty lattice, thus creating four bands in the first Bril- 
louin zone. For example let us consider a phase potential 
that varies between and 20o in n: 



20o, n : even 
0, n : odd 



(6) 



Note that this latter phase potential has the same 
strength as the one used before. In this latter case, by 
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using the same ansatz of Eq. (j3j) we directly obtain the 
dispersion relation corresponding to the new potential of 
Eq. ©. 



cos (2(Q + </>o)) = 8cos 2 (0 — 4>q) — Scos((f)o)cos(0 — 0o) 

+4cos 2 (0 o ) -3 



(7) 



A close examination of Eq. ([7]) reveals that this latter 
dispersion curve is identical to that of Eq. ((4]) , apart from 
a phase shift in both and Q. More specifically Q has 
shifted by an amount of — 0o while by </>o. Figure [5] 
shows a plot of this dispersion relation for 0o = 0.27T. 
The shift of origin compared to Fig. H^b) is evident in 
this figure. 
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FIG. 5. (Color online) Band structure of an optical mesh 
lattice with a non- symmetric step-like phase potential alter- 
nating between and 20o while 0o = 0.27T. Compared to 
the case of a symmetric phase potential (Fig. HJb)) the band 
structure is shifted from the center. 

In the rest of this work we consider for simplicity 
symmetric phase potentials for which the band structure 
is symmetric around Q = = 0. 



IV. OPTICAL DYNAMICS IN MESH LATTICES 

In this section we investigate optical dynamics in pas- 
sive mesh lattices. The impulse response of the system 
is of particular importance since is known to excite the 
entire band structure. For this reason only one of the 
waveguide elements is excited at m = 0. In what follows, 
the impulse response will be studied by using a[J = 1 with 
all the other elements in the array initially set to zero. 

Figure [6^a) shows the impulse response of this array 
lattice when </>o = 0. According to this figure light 
transport in this system exhibits upon spreading a 
highest slope of Q m ax = il/V^ with respect to the 
longitudinal axis. As we will see this result will be 
formally justified by considering the group velocity in 




FIG. 6. (Color online) Impulse response of a mesh lattice 
where the intensity profile of a™ and b™ is plotted; (a) 0o = 
(empty lattice), (b) 0o = OAtt. In both cases a® = 1 and all 
other elements are initially set to zero. 



this arrangement. The impulse response of the mesh 
lattice in the presence of a periodic phase potential 
with cj)Q = OAtt is also plotted in Fig. [6jb). In this last 
case, it becomes clearly apparent that the maximum 
speed of the excitation spreading becomes slower when 
0o increases. As in waveguide arrays [6], the impulse 
response can be viewed as a "ballistic" transport across 
the array. 

The band structure can also provide useful informa- 
tion concerning the evolution of more complicated initial 
excitations like localized wavepackets. More specifically, 
we consider initial distributions of and of the form 
f n e % Q° n where f n is a slowly varying envelope function 
(with a narrow spatial spectrum) and e^° n is a rapidly 
varying phase term signifying the central Bloch momen- 
tum Qo of this wavepacket. Therefore the propagation 
process of this discrete beam excitation can be effectively 
treated through a Fourier superposition of the Floquet- 
Bloch modes e ' l Q n e iern assumed before to analyze this 
system. In this regard, both the group velocity and the 
dispersion broadening of this wavepacket can be obtained 
by expanding the propagation constant 9 in a Taylor se- 
ries around Qo, that is: 



d# 



d 2 



= 0o + dQ^ {Q - Qo) + ^2 ho W - Qo) (8) 



dQ 2 



As in continuous lattices, the tangent of the beam angle 
(or "group velocity") is associated with the term: 



(9) 



d^ 

dQ lQo 

Using the dispersion Eq. (J2J), this group speed can then 
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be written as: 
dQ 



sin(2Q) 



4 [sin(20) - cos(<t>o)sin(0)] 



(10) 



where in this relation could be replaced from the dis- 
persion relation of Eq. ((5]) to obtain the right hand side 
as a function of Q and the band under consideration. Us- 
ing similar arguments, the discrete diffraction factor can 
be obtained from: 



D 



d 2 

d Q2 iQo 



(11) 



Figure depicts the beam angle £1 for several lattices 
with different amplitudes of the phase potential, </>o- Ac- 
cording to this figure, in an empty lattice (<fo = 0) 
this beam angle is zero at the center (Q = 0) and it 
is maximum at Q = = in the folded Brillouin zone 
scheme where to first order the dispersion relation dic- 
tates that Q = ±y/26. Hence, as previously indicated, 
the maximum slope expected in this configuration is 
Qmax — =L1/a/2. On the other hand for a lattice hav- 
ing a periodic phase potential, each band exhibits a zero 
group velocity at the center and at the edges (Q = ±7r/2) 
of the zone while the maximum happens somewhere in 
between. For the special case of 0o = 7r/2 the bands are 
translated in and hence in groups of two have identical 
group velocity curves, and as shown in Fig. |7fc) they lie 
on top of each other. 
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FIG. 7. (Color online) Beam tangent angle (f2) for several 
cases; (a) empty lattice (note that in this case the curve is 
folded to the reduced Brillouin zone), (b) for a lattice in the 
presence of periodic phase potential with 0o = 0.27T, (c) 0o = 

0.57T, (d) 00 = 0.77T. 



where 2A represents the Gaussian beamwidth and Qo 
designates the initial tilt in its phase front or central 
Bloch momentum. In this case the same input profile 
is assumed for 6^ in order to symmetrize the dynamics. 
Figure [8] shows the propagation dynamics of this Gaus- 
sian beam in this mesh lattice. Here the lattice involves a 
periodic phase potential with Qo = 0.27T. The Gaussian 
beam width 2A is large enough to avoid the diffraction 
effects and in addition its tilt is Qo = 0.257T. Accord- 
ing to this figure four independent beams (of the same 
width) result from this initial excitation, each emanating 
from a corresponding band, and propagating in different 
directions. To elucidate these results, the band structure 
is also plotted in this same Fig. [8jc) where the arrows 
perpendicular to the bands indicate the propagation di- 
rection of each of these four beams. 
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FIG. 8. (Color online) Gaussian wavepacket propagating in 
a mesh lattice. The beam has a width of 2A = 30 and an 
initial phase tilt of Qo = 0.257T. The lattice has a phase 
potential of 0o = 0.2tt (a) intensity |aC| 2 , (b) intensity of 
|&tT| 2 j ( c ) band structure of the lattice with the dashed line 
crossing the band at four points at Qo = 0.257T and the arrows 
show the propagation direction of the four resulting beams, 
(d) intensity profile of the initial Gaussian beam, (e) \aif\ 
intensity profile of a™ at the last discrete longitudinal step 
(here M = 300), (e) \b?f | 2 . 



To demonstrate some of these transport effects, let 
us consider for example the evolution of a Gaussian 
wavepacket having the following initial profile: 



a = e -(n/A) 2 e iQ n 



(12) 



Finally in order to investigate diffraction effects in pas- 
sive mesh systems, we consider the propagation proper- 
ties of a relatively narrow Gaussian wavepacket. Figure [9] 
depicts the propagation dynamics of a Gaussian beam 
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with a width of 2A = 8 in a lattice with 0o = 0.57T. The 
figures compare the beam propagation for two different 
values of Qo, and 0.257T. According to this figure when 
Qo = 0, the beam has a very low transverse velocity 
and experiences a considerable degree of diffraction. As 
shown in the other panels, when the beam is launched at 
the dispersion free point of the band (Qo = 0.257r), where 
D = and the transverse group velocity is maximum, the 
diffraction effects are negligible. 




FIG. 9. (Color online) Diffraction properties of a Gaussian 
beam in a mesh lattice with 0o = 0.57T. The Gaussian beam 
has a width of 2A = 8 while the initial phase tilt is: (a) 
Qo = 0, (b) Qo = 0.25tt. 

According to Fig. |4] this selection of 0o leads to four 
bands. Figure [9ja) depicts Gaussian beam spreading 
at Qo = and at the same time interference effects 
resulting from the excitation of multiple bands. On 
the other hand for Qo = 0.257T two Gaussian beams 
symmetrically emerge with two different propagation 
speeds. Yet, the interference pattern in each of the two 
branches demonstrates that all four bands are actually 
in play in these dynamics. Notice however that at 
this point little beam spreading occur since for these 
parameters D = 0. 



V. PT- SYMMETRIC BUILDING BLOCK 

Before exploring a large-scale PT-symmetric mesh lat- 
tice, it is worth analyzing the elemental building block 
involved in such a network. Figure [TUTa) shows a VT- 
symmetric coupler where the gain and loss is uniformly 
distributed along the two arms, a structure similar to 
that considered in previous experimental studies [22|, [23| . 
Figure [TDTb) on the other hand depicts a passive coupler 



where the gain and loss mechanisms are separately in- 
serted in the two arms only. Here we show that this new 
type of T^T-symmetric coupler displays exactly the same 
behavior and characteristics of a standard PT-coupler 
arrangement considered before. 




FIG. 10. (Color online) A distributed T^T-symmetric coupler 
and a "PT-synthetic coupler; (a) The "PT-coupler is composed 
of two similar dielectric waveguides coupled to each other, 
with one experiencing gain (red) while the other an equal 
amount of loss (blue), (b) A PT- synthetic coupler is com- 
posed of a passive coupler while the gain and loss waveguides 
are separately used in the arms. 

In Figure [TOT b) we assume a 50:50 passive directional 
coupler connected to two arms, one providing amplifica- 
tion (red) while the other an equal amount of loss (blue). 
We assume that each arm delivers an amplification or 
attenuation of e ±7 / 2 right before and after the coupler. 
Hence the modal amplitudes a' and b' at the output of 
this system, are related to those at the input ports, a and 
6, in the following way: 

a'\_(e + ^ 2 \J_/1 i\ (e +l/2 \ (a\ 
b>)-\ e-^J V2 V V V e-^V \b) 

(13) 

in which case 

where a and b represent optical amplitudes in the gain 
and loss channels respectively. The two supermodes and 
their respective eigenvalues of this system can be readily 
found. Depending on the amount of gain/loss in the sys- 
tem two regimes can be distinguished; if 7 < cosh -1 (\/2) 
this VT system is operating below the PT-symmetry 
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breaking threshold and its supermodes are given by: 



± e ±iri 



(15) 



where cos(uo) = -j=cosh(j) and sin(uo) = -^=cos(rj). 

Thus for 7 < cosh~ 1 (^/2) the two modes repeat them- 
selves after passing through this discrete system except 
from a trivial phase shift of ±cj. On the other hand 
if 7 > cosh -1 (V2) the system operates above the VT- 
symmetry breaking threshold and: 



(16) 



where cosh(uo) = -^=cosh(j) and sinh(uo) = -^=sinh(rj). 
Interestingly this same behavior is displayed by a stan- 
dard T^T-symmetric coupler where the gain and loss is 
continuously distributed. Finally at exactly the VT- 
symmetry breaking threshold 7 = cosh~ 1 (^) the two 
supermodes collapse to one and thus: 



(17) 



which clearly shows the existence of a phase difference 
of 7r/2 between the two waveguides. 

It is worth noting that this arrangement has certain 
advantages over a standard distributed PT-symmetric 
coupler. First of all it is experimentally easier to achieve 
the delicate balance required for VT symmetry. In 
addition the coupling and amplification/attenuation 
process take place in two separate steps so there are no 
physical restrictions imposed by the Kramers-Kronig 
relations. As previously mentioned, these effects have so 
far hindered progress in implementing large-scale VT- 
symmetric networks, since they limit the possibility of 
achieving the required values for gain/loss and refractive 
index, all at the same time. 



VI. PT- SYNTHETIC MESH NETWORKS 

Figure ITTf a) shows a VT- symmetric mesh lattice made 
of PT-synthetic couplers, identical to that of Fig. fTQlb). 
In addition phase elements are inserted in this same 
lattice (shown by circles in Fig HTT a)) in order to pro- 
vide the needed real part in the potential function. In 
Fig. fTTTT b) the distributions of phase modulation and 
that of gain/loss are plotted as a function of the discrete 
position n - clearly satisfying the requirement for VT- 
symmetry, i.e. an even distribution for the phase and an 
odd distribution for the gain/loss profile in n. In fact 
a comparison with continuous systems suggests that the 
phase and gain/loss in discrete elements play the role of 
the real and imaginary parts in the refractive index re- 
spectively. By considering an amplification/attenuation 



factor of e ±7 / 2 in each waveguide section between two 
subsequent couplers, then one can show that light prop- 
agation in this PT-synthetic mesh network is governed 
by the following discrete evolution equations: 



(a) 



(b) 
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FIG. 11. (Color online) (a) A 7^ T- synthetic mesh lattice, (b) 
transverse distribution of the phase potential (symmetric) and 
gain/loss (antisymmetric). 
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To understand the behavior of this system, the band 
structure should be first determined. By adopting the 
same ansatz of Eq|3j one can derive the following disper- 
sion relation for this VT lattice: 



cos(2Q) = 8cos 2 (6) - 8cosh(j)cos((j)o)cos(0) 
+4cos 2 (0o) - 4 + cosh(2j) 



(19) 



Figure [12] shows the band structure of this system for 
several different values of the phase potential amplitude 
0o and gain/loss coefficients 7. In each case the real parts 
of the propagation constant (9) is plotted in blue while 
the imaginary parts are shown in red. 

As it is illustrated in this figure, the presence of a 
symmetric phase potential in this system tends to pull 
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FIG. 12. (Color online) Band structure of VT synthetic mesh 
lattice for several values of 0o and 7. In these plots the real 
part of propagation constant, is indicated in blue, while the 
imaginary part in red. 



by: 

7 = cosh' 1 (2cos((j)o) - Vcos(20 o )) (21) 

This relation dictates the merging condition for the first 
two bands and is only valid for < </>o < 7r/4, consistent 
with our previous observations. To find the correspond- 
ing relation for the band merging occurring at the edges, 
in Eq. (fT9|) we set Q = tt/2 , which in turn leads to a sec- 
ond order algebraic equation in cos(6). Since we expect 
that the two eigenvalues will collapse into one (excep- 
tional point), one may use this degeneracy condition in 
Eq. (fT9|) at Q = tt/2. After setting the discriminant of 
the quadratic equation to zero one finds that: 

7 = cosh' 1 (Vfj « 0.8814 (22) 

This last relation provides the PT-threshold for band 
merging at the edges of the Brillouin zone and is indepen- 
dent of 0o- Interestingly this same value 7 = cosh' 1 (\/2) 
coincides with the critical PT-threshold of the basic 
unit involved in this lattice, as found in section V. 



apart the bands thus creating a band gap, while the 
antisymmetric gain/loss tends instead to close the gap. 
The system is said to be operating below the VT- 
symmetry breaking threshold as long as the eigenvalues 
associated with all bands are real. However at a critical 
amount of gain/loss the bands merge at the so called 
exceptional points, and for even higher gain/loss values, 
sections with conjugate imaginary eigenvalues appear in 
the bands. 



Figure [13] depicts the PT-symmetry breaking thresh- 
old in the parameter space of </>o and 7. The area below 
the curve corresponds to the case where the system oper- 
ates in the exact VT phase where all the eigenvalues are 
real. On the curve symmetry breaking occurs and above 
this line the spectrum is in general complex. The top 
fiat line of this curve corresponds to the critical value of 
0.8814 while the part between < </>o < 7r/4 can be ob- 
tained from Eq. ([2T]) . The other segement symmetrically 
follows. 



In what follows, we consider the case where </>o is 
fixed and discuss how the band structure will change by 
gradually increasing the gain/loss coefficient 7. Analysis 
shows, that for a given value of 0o, the first band merg- 
ing occurs at two different positions; if < </>o < 7r/4, the 
bands merge at Q = 6 = and the second band gap re- 
mains open till reaching a critical value of gain/loss coeffi- 
cient 7. For even higher gain/loss values the system finds 
itself in the broken phase regime. For 7r/4 < 0o < 
on the other hand all bands are open till a critical point. 
Exactly at this threshold, the band gap at the edges 
of the Brillouin zone at Q = ±7r/2 closes while the 
first band gap remains open till reaching another criti- 
cal point where it eventually evaporates. Based on these 
observations analytical results for the symmetry breaking 
point can be obtained. We first consider the case where 
< </>o < 7r/4. In this case, as 7 increases, we expect 
that for a fixed 0o, the symmetry breaking will occur at 
Q = 6 = 0. Therefore Eq. ([T9]) can be rewritten as: 
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7T/4 tt/2 3n/4 n 

( Po 

FIG. 13. (Color online) PT-symmetry breaking threshold 
curve in a two dimensional parameter space of 0o and 7. The 
region below the curve corresponds to the exact PT-phase 
while the region above the curve designates the domain where 
VT symmetry is broken. 



cosh 2 (~/) - 4cos((f)o)cosh(j) + 2co5 2 (^ ) + 1=0 (20) 



From here one can easily show that this critical 7 is given 



To dynamically explore the symmetry breaking thresh- 
old, the impulse response of our system is studied. Since 
the impulse is expected to excite the entire band of 
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this mesh lattice, one should expect that an exponen- 
tial growth in the total energy of the system should be 
observed once the PT-symmetry is broken. 
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FIG. 14. (Color online) Impulse response of the PT- 
symmetric lattice with a periodic phase potential of 0o = 0.2tt 
while several different amounts of gain/loss are considered; (a) 
7 = (the passive lattice), (b) 7 = 0.3 (below threshold), (c) 
7 = 0.35 (at threshold), (d) 7 = 0.4 (above threshold) 

Figure HH shows the impulse response (a° = 1 , while 
all other elements are initially zero) of a PT-symmetric 
mesh lattice for several different values of gain/loss 7 
when <j)o = 0.2tt. This range covers the passive scenario, 
or the case where the system operates below, at, and 
above the PT-symmetry breaking threshold. The total 
energy in the system E m = |a™| 2 + |6™| 2 , is also 
plotted in each case at each discrete step of propagation, 
m in Fig. [T4J While for the passive system (7 = 0) the 
total energy remains constant during propagation, for 
a PT-symmetric lattice used below its threshold the 
total energy tends to oscillate during propagation -but 
always remains below a certain bound. Note that such 
power oscillations were previously encountered in other 
'PT-symmetric periodic structures [l8j|. At exactly the 
PT-threshold a linear growth in energy is observed-see 
Fig. fT4f c). Finally above threshold an exponential 
growth in energy is observed as expected from a system 
involving complex eigenvalues- Fig. fT4f d). 

To further explore the behavior of this PT-synthetic 
mesh lattice, we use at the input a Gaussian wavepacket, 
as in Eq. ([T2]) . Indeed by exciting this system with a 



wide input beam (that has a narrow spectrum) one can 
selectively excite different sections of the band structure. 
We now consider a PT-symmetric mesh lattice with a 
periodic phase potential of amplitude c^o = 0.2tt and a 
gain/loss factor of 7 = 0.4. The band structure corre- 
sponding to this structure is plotted in Fig. [T2l 
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FIG. 15. (Color online) Gaussian beam propagation in a PT- 
symmetric lattice operating in the broken PT phase regime. 
The lattice has a periodic phase potential of amplitude 0o = 
0.27T and a gain/loss factor of 7 = 0.4. The Gaussian beam 
has a width of 2A = 30 and is launched with three different 
values of initial phase tilt; (a) Qo = 0, (b) Qo = 0.257T, (c) 
Qo = 0.57T. In (a) the intensities are only shown up to a level 
of 100. 

Figure [15] depicts the propagation of a Gaussian 
wavepacket in this lattice, when launched with a Bloch 
momentum Qo- Three different values for Qo have 
been selected for this example: Qo = 0, 0.257T and 
0.57T. According to Fig. [15] while for the first case an 
exponential energy growth is observed, for the other two 
cases energy remains essentially bounded. These results 
reveal that even above the PT-symmetry breaking 
threshold, non-growing/decaying modes can be excited 
in such systems. This all depends on which section of 
the band structure is excited by the initial conditions. 

Compared to a passive mesh lattice, the band struc- 
ture of its PT-symmetric counterpart reveals another in- 
teresting property. As previously discussed, the maxi- 
mum beam transport angle (flmax) in an empty lattice 



11 



is 1/ a/2, and even in the presence of a periodic phase po- 
tential this angle is always less than this maximum trans- 
verse velocity. However according to the Fig. [I2j when 
approaching the exceptional points from the real section 
(blue part) of the band, its slope tends to considerably 
increase and eventually approaches infinity around the 
exceptional points. 
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FIG. 16. (Color online) A broad Gaussian beam propagat- 
ing in a passive and a PT-symmetric lattice; (a) evolution 
of the Gaussian beam in a passive empty lattice, (b) in a 
PT-symmetric lattice, (c) normalized intensity profiles of the 
beam at the last propagation step (m — 300) in both lattices. 
The parameters of the PT lattice are 7 = 0.039 and 0o = 0. 
The Gaussian beam has a beam width of 2A = 400 and an 
initial phase front tilt of Qo = 0.98177T. 

Figure [16] compares the propagation of a Gaussian 
beam in a passive and a PT-symmetric mesh lattice 
operating above threshold. Both lattices are excited with 
the same Gaussian beam having a Bloch momentum 
Qo? which is chosen to be close to the exceptional point 
of the PT-symmetric lattice. Close to this exceptional 
point, the slope of the band structure tends to infinity 
and therefore, the associated group velocity can become 
almost arbitrarily high for any narrow-bandwidth 
wavepacket. While the maximum beam angle in passive 
empty lattice is ~ 0.7 (which is close to the maximum) 
for the PT-symmetric lattice this angle is approximately 
1.04 which is certainly above the maximum limit of the 
passive lattice. This effect has in fact a counterpart in 
continuous media. As previously shown, in the presence 
of a gain medium [54], [55j and in PT-symmetric and 
gain/loss gratings and lattices [38|, [EH, use d close 
to the exceptional points, the group velocity of light 



can be superluminal. It should be noted however 
that none of these effects violates causality since non 
causal waveforms are used for excitation. Indeed, this 
superluminal propagation of the intensity peak is en- 
abled by a gain-assisted growth of the distribution's tails. 

Finally we investigate the concept of unidirectional in- 
visibility in a PT-symmetric mesh lattice. As recently 
predicted [35|, [36[ PT-symmetric periodic structures like 
gratings can exhibit surprising behavior like unidirec- 
tional invisibility and intriguing reflection characteristics. 
More specifically, light propagating in such a system can 
experience reduced or enhanced reflections depending on 
the direction of propagation. Even more remarkable is 
what happens right at PT threshold: in this case light 
waves entering this arrangement from one side do not ex- 
perience any reflection and can fully traverse the grating 
with unity transmission. Given that this occurs without 
acquiring any phase imprint from this PT system, the 
periodic structure is essentially invisible. Similar effects 
also occur in PT-symmetric mesh lattices. This can be 
demonstrated in a system where one lattice is embedded 
in another lattice. This is achieved using for example the 
following phase modulation potential: 



0, |n| > n 
+</>o, \n\ < no, n : even 
—0o, \n\ < no, n : odd 



(23) 



where 2no is the width of this PT grating lattice. In 
this same way the anti-symmetric gain/loss profile is 
also imposed only in the grating region |n| < no- 

To investigate this latter process, a PT-symmetric lat- 
tice having 2no = 40 layers is embedded inside an empty 
lattice. The empty lattice is excited with a Gaussian 
beam (as in Eq. ([12]) ) with an initial phase tilt of Qo = it. 
Scattering of the beam from the left and right side of this 
PT-symmetry grating is depicted in Fig. fT71 

Figure [TTl shows the scattering of Gaussian wavepacket 
from this grating when 2no = 40, 0o = 0.27T, and 
7 = 0.3507. The extent of this grating is shown by 
the dashed lines. Figure fT7t a) shows how the Gaussian 
beam is scattered or transmitted by this grating when 
7 = 0. The effects of reflection and reduced transmission 
are evident in this figure. In Fig. [TTIf b) we show these 
same dynamics for 7 = 0.3507 and when the Gaussian 
excites the grating from the left. Both the a™, b™ 
channels are excited during this process. In this case 
no reflections occur when the grating is used close to 
the exceptional point. Essentially in this regime the 
grating leaves no mark on the beam itself and hence is 
practically invisible. Note that the associated splinter 
beams in this figure do not represent reflections-they 
simply emerge from the two different bands associated 
with the empty lattice. On the other hand, Fig. fTTT c) 
shows what will happen when the Gaussian beam excites 
the right side of this same PT-symmetric lattice grating. 
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FIG. 17. (Color online) Unidirectional scattering of a Gaus- 
sian beam from a PT-symmetric grating. A PT-symmetric 
grating having 40 layers with a phase potential of </>o = 0.2tv 
and a gain/loss of 7 = 0.3507 (corresponding to VT thresh- 
old) is established in the middle of an empty lattice. The 
dashed line shows the grating region, (a) Scattering from 
the passive grating, (b) left side scattering from this VT- 
symmetric grating, (c) right side scattering from this same 
grating. 



In this latter case pronounced reflections occur (even 
exceeding unity) and the grating ceases to be invisible 
to light. 



amplifiers/attenuators can be employed to realize the 
antisymmetric imaginary part of the VT potential. 
What makes these optical lattices different from previ- 
ously known waveguide array versions is the presence 
of discreetness in both the transverse and longitudinal 
directions. The band structure of these systems has 
been systematically analyzed and an analytic expression 
was obtained for their dispersion relation. We have 
shown that the band structure is periodic in both the 
propagation constant and transverse Bloch momentum 
while the Brillouin zone of these lattices displays in 
general four bands. In addition we found that the 
shape of the bands and band gaps can be effectively 
controlled using phase elements. Interestingly, through 
a proper phase modulation the band structure can be 
arbitrary shifted from the center. It should be noted that 
shifting the bands in standard optical waveguide array 
systems is not straightforward and typically requires the 
presence of external magnetic field effects. The impulse 
light dynamics as well as wavepacket excitations were 
numerically explored and related to the properties of 
the band structure. The elementary PT-building block 
involved in these mesh arrangements was examined 
and its symmetry breaking threshold was determined. 
As it was discussed, what could greatly facilitate the 
physical realization of such a large scale PT-symmetric 
mesh lattice is the fact that couplings and amplifica- 
tion/attenuation can be independently controlled within 
the basic building block of this lattice. Band merging 
effects in these lattices were investigated and the con- 
ditions for spontaneous PT-symmetry breaking were 
explicitly obtained in terms of relevant parameters. The 
response of these systems under impulse and broad beam 
excitation was investigated in terms of their respective 
band structure. As it was shown, light dynamics in 
the PT-symmetric lattice exhibits certain peculiarities 
that are otherwise impossible in its passive counter- 
part. These include for example power oscillations 
and transitions from neutral to exponentially grow- 
ing regimes. The possibility of superluminal transport 
along with unidirectional invisibility was also considered. 



VII. CONCLUSIONS 

In this work we have studied the properties of a new 
class of periodic structures both in the passive as well in 
the PT-symmetric regime. These optical mesh lattices 
are in essence waveguide arrays that are discretely and 
periodically coupled to each other along the propagation 
direction. In addition phase elements can also be used 
in appropriate positions to control the phases while 
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